Self- Accelerating Universe in Galileon Cosmology 
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We present a cosmological model with a solution that self-accelerates at late-times without signs 
of ghost instabilities on small scales. The model is a natural extension of the Brans-Dicke (BD) 
theory including a non-linear derivative interaction, which appears in a theory with the Galilean 
shift symmetry. The existence of the self-accelerating universe requires a negative BD parameter 
but, thanks to the non-linear term, small fluctuations around the solution are stable on small scales. 
General relativity is recovered at early times and on small scales by this non-linear interaction via the 
Vainshtein mechanism. At late time, gravity is strongly modified and the background cosmology 
shows a phantom-like behaviour and the growth rate of structure formation is enhanced. Thus 
this model leaves distinct signatures in cosmological observations and it can be distinguished from 
standard ACDM cosmology. 



I. Introduction 

Observational evidence indicates that the universe is 
undergoing a period of accelerated expansion [l|, 0, d, 0] • 
This late-time acceleration is one of the biggest open is- 
sues in cosmology today. This would mean that 70% 
of the energy content of the universe is dominated by 
unknown dark energy. The best candidate is vacuum en- 
ergy, but it is more than 50 orders of magnitude larger 
than the required value of the cosmological constant to 
explain the current acceleration. Alternative explana- 
tions are therefore actively being pursued. 

Modified Gravity is an alternative that interprets this 
acceleration as a sign that our knowledge of how gravity 
works at large distances is incomplete and tries to ex- 
plain the acceleration by a modification of General Rela- 
tivity (GR) . One of the most studied attemptSjthe Dvali- 
Gabadadze-Porrati (DGP) braneworld model [5, 6], has a 
self- accelerating solution, i.e. a solution in which cosmic 
acceleration arises even when no matter is present in our 
universe Unfortunately this solution was found to 
be plagued by ghost instabilities that cast doubt on the 
validity of the model (see [8| for a review and references 
therein). Other similar alternatives were found (see e.g. 
@, 03) but it seems incredibly difficult to separate the 
presence of the ghost from the acceleration [ll| ■ 

Recently, an infrared modification of gravity was pro- 
posed, which is a generalization of the 4D effective theory 
in the DGP model [l2j (see also [IH). A novel feature 
of this effective theory is that the theory is invariant un- 



der the Galilean shift symmetry, the constant shift of the 
gradient of the scalar field — > <9 M </> + c^. which keeps 
the equation of motion at second order [14| . Then they 
found self-accelerating de Sitter solutions with no ghost- 
like instabilities. However, their analysis is valid only 
for weak gravity in flat space-time and this result might 
change when the model is covariantized [IH m, m. in 
fact it was shown that the Galilean symmetry cannot be 
preserved once the theory is covariantized but it is still 
possible to keep its desired properties. For example, the 
equation of motion for the scalar field can remain of sec- 
ond order, which is essential because the higher derivative 
theories contain extra degrees of freedom that are usually 
related to instabilities. 

In this letter, we demonstrate that it is possible to con- 
struct a covariant model that keeps desired properties of 
the Galileon model, i.e. the existence of self-accelerating 
universe with no ghost-like instabilities on small scales. 
We consider the following action (see also 



<t>R - -(D4>) 2 + c4 + /(0)D0(D0) 2 + C n 

<P 

(1) 

where (Dtp) 2 = D a 4>D a (j> and C m is the matter la- 
grangian. The cubic interaction is the unique form of 
interactions at cubic order that keeps the field equation 
for 4> of second-order pjj]. The Einstein and field equa- 
tions are given by 
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where T^ v is the energy-momentum tensor for matter. 
II. Background Cosmology 

We found that the simplest choice of / to obtain the 
self- accelerated universe is given by /(</>) = I/A/ 2 2 , 



where M is the parameter of the model. We can also 
redefine c as c — — 2A so that it acts like a cosmological 
constant. For Friedmann-Robertson- Walker spacetime, 
the Einstein and field equations, and ^ give 
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where H is the Hubble parameter, p is the energy density, 
p is the pressure and P = 4>/4>. 

Self- Accelerating solutions: A self-accelerating solution 
can be found by putting A = p = p = and searching for 
a solution with H = P = 0. From Eq. ([5|) one can easily 
show that this solution must satisfy 



(7) 



For this solution to exist, the BD parameter must sat- 
isfy ui < —4/3. The Friedmann and field equations now 
yield 
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This is a self-accelerating dc Sitter solution in which 
the linear terms are equated by the non-linear terms 
coming from the cubic interaction in Eq. ^ . Given that 
H 2 has to be positive, we have to pick the negative sign 
branch in Eq. (J7J). In order to describe the acceleration 
today, M should be fine-tuned M ~ Hq where H$ is the 
present-day Hubble parameter. 

Early-time solutions: To recover GR at early times P 
should be smaller than H. In fact, GR is recovered 
if P = 0. This is realized at high energies when the 
non-linear terms in the field equation dominate over 
the linear term. In this case, from Eqs. (JUIHI), we get 
P = M/V3 ~ Ho < H and we indeed recover the GR 
solution. This is a cosmological version of the Vainshtein 
mechanism to screen the scalar field [l6j . 
Numerical solutions: We can solve Eqs. (JSJIS]) numeri- 
cally with appropriate initial conditions. Throughout 
the paper, a flat cosmology is assumed but the in- 
clusion of spatial curvature is straightforward. We 
set initial conditions at matter domination era and 
found the appropriate set of initial conditions in or- 
der to get acceleration at present time (defined by 



^m(io) = p(to)/3H(t ) 2 (f)(to) — 0.3 as an example). 
Note that Newton's constant is given by 8ttG — l/4>{ti) 
in this model where ti is the initial time in the matter 
dominated era. We obtained cosmological solutions that 
mimic ACDM, i.e. the scalar field becomes important at 
late-times and the universe enters a period of accelerated 
expansion, approaching the self-accelerating solution. 
We also observed that whenever P(U) is appropriately 
small, P(U) <C H(ti) the result is insensitive to the ini- 
tial value of P(U). In order to have f2 m (to) = 0.3 today, 
the parameter M must be fine tuned as M = O.I5i?o 
for lu = -50 and M = 0.024iJ for lo = -500. The 
contributions of P in Eq. ^ can be combined into an 
effective dark energy, such that the Friedmann equation 
takes the form H 2 = (8irG) (p + p c ff) /3. This allows 
us to plot the effective equation of state, u> c ff, of this 
effective dark energy for different values of the BD 
parameter as a function of redshift (Fig. []}. We also 




FIG. 1: The effective equation of state iu e a as a function of 
redshift for various values of the Brans-Dicke parameter. 

computed the comoving distance r(z) = J z dz 1 / H(z'), 
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which is plotted in Fig. [2] Due the phantom-like be- 
haviour w e g < — 1, the distance is larger than in ACDM. 
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FIG. 2: The comoving distance r(z) as a function of redshift 
for different values of to. 



III. Perturbations 

Stability: We now study the stability of the model against 
small perturbations (f> — > 0(1 + ip). To do this we ex- 
pand the field equation ([6]) at linear order. We consider 
perturbations on small scales where we can neglect the 
expansion of the universe and the effect of metric pertur- 
bations. The evolution equation for tp becomes 



d t {t)(p + d x {t)V 2 p = 0, 



where 



dt{t) = 3 + 2cj 
d x (t) = 3 + 2lu+^ 
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In the usual BD theory, if the BD parameter is smaller 
than —3/2, then d± is negative and ip becomes a ghost. 
In our model, the non-linear interaction term changes 
the sign of d t - We found that as long as P is positive, 
dt and d x are positive at all times. This can be realized 
if u> < —2 (see Eq. (0). One can also calculate the 
sound speed of these perturbations, c 2 = d x /d t - At early 
times, the terms proportional to H dominate and thus 
c 2 s = 2/3, which is subluminal. The behaviour of c s (t) for 
the background solutions discussed above can be found 
in Fig. [31 The fluctuations can be superluminal for large 
| a; | at the transition from the matter dominated era to 
the accelerating phase and this would impose an upper 
bound for \uj\. 

Growth Factor. The evolution equation for the cold dark 
matter over-density 5 is given by: 



5 + 2H5 



(12) 
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FIG. 3: Plot of c a as a function of redshift for various val- 
ues of the Brans-Dicke parameter. As can be seen during 
the transition to the accelerating phase the perturbations go 
superluminal for u) < —190. 



where the perturbed line element is given by ds 2 = — (1 + 
2^)dt 2 + a(t) 2 (l + 2^)S ij dx i dx j . On small scales, the 
Einstein equations (|2|) give 
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These can be combined with the field equation (from 
©) 



V 2 
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where we use the quasi-static approximations and ne- 
glected the time derivative term. Then the equation for 
5 is obtained as 



S + 2HS = 4nG cB pS, 
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(17) 



The effective Newton's constant G c s is close to G at 
early times but it becomes larger at late times. This is 
because the scalar mode <p gives an additional attractive 
force. We should note that this is opposite to the DGP 
self- accelerating solutions where the effective Newton 
constant is smaller due to the fact that the scalar mode 
is a ghost and it mediates a repulsive force. Solving 
the evolution equation numerically we obtained the 
growth factor 8/a as is shown in Fig. |U Due to the 
enhancement of Newton's constant, the growth rate is 
enhanced compared with the ACDM case. 

IV. Discussion 

By including a cubic order derivative interaction in the 
action for the Brans-Dicke theory, inspired by the DGP 
effective action and the Galileon model [121. we found 
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model too, but this only occurs for large enough values of 
the parameter |a>| in our model. The full study of cosmo- 
logical perturbations is necessary to ensure the stability 
of the model on horizon scales and this is an important 
open issue. 

At late times, gravity is strongly modified. This raises 
a question whether it is possible to evade the strong con- 
straints on the deviation from GR such as the solar sys- 
tem tests. Again the cubic interaction term comes to 
rescue. On small scales, including the non-linear term, 
the equation for (p looks like 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

FIG. 4: The growth rate ^ as a function of redshift for dif- 
ferent values of the BD parameter, ui. 

a self-accelerating solution with no signs of ghost insta- 
bilities on small scales. This solution exists if the BD 
parameter is smaller than —2. Whenever the condition 
H ^> <j)/<f> > is satisfied initially, the resultant solutions 
are insensitive to the initial conditions for the scalar field. 
In order to get the observed acceleration, the parameter 
M needs to be fine-tuned, M ~ Hq. A remarkable prop- 
erty of this cubic interaction is that it helps to recover 
GR at high energies. At high energies, this cubic interac- 
tion term ensures </>/</> ~ M <C H. This is a cosmological 
verision of the so-called Vainshtein mechanism [l6| 

The kinetic terms for small fluctuations on small scales, 
dt and d x are positive indicating that the model is free of 
ghost and tachyonic instabilities. The cubic interaction 
again plays a crucial role here. Without the cubic inter- 
action, the fluctuations would be a ghost if to < —4/3, 
but the cubic interaction gives an additional contribu- 
tion to the kinetic term for the fluctuations that changes 
the sign of the kinetic terms. The sound speed can ex- 
ceed one depending on model parameters and this super- 
luminality would impose constraints on the model param- 
eters to and M. We should note that there is the issue 
of super-luminal propagations in the DGP and Galileon 



where tp <C 1. This is the same equation as the one for the 
brane bending mode in the DGP model (see [HI]). For 
a static spherically symmetric source, it has been shown 
that GR is recovered on small scales r < r* where the 
Vainshtein radius is given by ~ r g /M 2 d'? cl where 
r g is the Schwarzschild radius of the source. Since M~ l 
is larger than the cosmological horizon scales, the Vain- 
shtein radius is sufficiently large to recover GR and evade 
solar system constraints, but a small deviation from GR 
could be tested in future experiments. 

In summary the inclusion of the cubic interaction in the 
BD theory provides surprisingly rich phenomenology. It 
is then important to study the effects of the other higher 
order interactions discussed in [l2l |. Moreover, since the 
Galilean symmetry is lost in the covariant version of the 
theory, there is no way to prevent higher order interac- 
tions from being generated and it is crucial to find a way 
to control these terms. This is related to the strong cou- 
pling problem, which constitutes yet another problem in 
the DGP model. This remains an open issue. 
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